Off-lattice Noise Reduced Diffusion-limited Aggregation in Three Dimensions 
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Using off-lattice noise reduction it is possible to estimate the asymptotic properties of diffusion- 
limited aggregation clusters grown in three dimensions with greater accuracy than would otherwise 
be possible. The fractal dimension of these aggregates is found to be 2.50 ± 0.01, in agreement with 
earlier studies, and the asymptotic value of the relative penetration depth is — = 0.122 ± 0.002. 
The multipole powers of the growth measure also exhibit universal asymptotes. The fixed point noise 
reduction is estimated to be e ~ 0.0035 meaning that large clusters can be identified with a low 
noise regime. The slowest correction to scaling exponents are measured for a number of properties of 
the clusters, and the exponent for the relative penetration depth and quadrupole moment are found 
to be significantly different from each other. The relative penetration depth exhibits the slowest 
correction to scaling of all quantities, which is consistent with a theoretical result derived in two 
dimensions. 



CO 

03 



i 

C 

O 
o 



> 

O 
in 

O 

o 

C3 



I 

-a 
c 

o 
o 



X 



Diffusion-limited aggregation (DLA) |l| is an exten- 
sively studied model of diffusion limited growth which 
appears to capture the essential features of many differ- 
ent physical growth phenomena 0,0, El HI- However, the 
fractals generated have evaded a complete understanding 
for many years and there has recently been controversy 
over whether diffusion-limited aggregates are truly frac- 
tal @. 

DLA is modelled p) by allowing particles to randomly 
walk from a sphere surrounding the cluster, one at a time, 
until they contact the cluster, at which point they are ir- 
reversibly stuck. Detailed study |j| has shown that DLA 
growth in two dimensions does approach true fractal scal- 
ing, but with slowly decaying corrections to scaling of the 
form 



Qn = Qc 



CN~ 



(1) 



where Qn is some property of the cluster, tending to- 
wards the value Qoo as the number of particles in the 
cluster, N, tends to infinity. Here the correction to scal- 
ing exponent v is expected to exhibit some universality 
whilst the constant C will not. For aggregates grown in 
two dimensions, it has been argued Q that there should 
be no quantity whose correction to scaling is slower than 
that of the relative penetration depth, „ , where Rdep 
is the average radius at which new particles are deposited 
and £ is the standard deviation of the same. Here and 
below we take as origin the centroid of the depositing 
particles. 

When studying DLA grown on a lattice, reducing the 
shot noise associated with the growth being by discrete 
particles Q has proved valuable in understanding the 
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asymptotic properties of clusters, complicated by their 
sensitivity to lattice anisotropy [3, Using a con- 

formal map from the unit circle to the boundary of a 
growing cluster, Hastings & Levitov introduced a 
technique for implementing a noise reduction scheme for 
DLA clusters grown off-lattice. Rather than adding par- 
ticles to the cluster, bumps were added to the conformal 
map. Ball et al. [H built on this work, allowing crescent- 
shaped bumps to be added to DLA clusters without the 
need for a conformal map. In this approach, the par- 
ticles are allowed to diffuse normally until they contact 
the cluster. At this point the new particle is touching the 
cluster, and the distance between the centre of the new 
particle and the centre of the particle it contacted in the 
cluster is equal to the diameter of one particle. To imple- 
ment noise reduction, this distance is reduced by a factor 
A < 1, so that the new particle is deposited partially in- 
side the cluster: the effect is to protrude a shallow bump 
of height A on the cluster perimeter. Since this method 
does not rely on a conformal map it allows the growth of 
noise reduced DLA clusters off-lattice in any dimension. 

Most of the work on DLA has been restricted to two 
dimensions. Meakin pioneered work on DLA in higher 
dimensions, growing clusters in dimensions up to d = 8 
[l3j. Much of that work has focused on estimating the 
fractal dimension of DLA clusters and the scaling of the 
relative penetration depth 0, 0, , yet firm conclu- 
sions have proved difficult. There has also been some 
progress measuring the multifractal spectrum of DLA in 
three dimensions |I^lia . ll9| . However, Davidovitch et al. 
|2C| recently claimed that all previous attempts to mea- 
sure f(a) in two dimensions are poorly converged, so the 
early three dimensional measurements should be taken 
with caution. Other work has also examined DLA on a 
cubic lattice in the limit of zero noise j^J, and the ex- 
tension of the fixed scale transformation to 3 dimensions 

In this paper we exploit the new noise reduction tech- 
niques to explore the convergence to scaling of DLA in 
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three dimensions. We grew 1000 DLA clusters off-lattice 
in three dimensions spanning five different values of the 
noise reduction A = 1, 0.3, 0.1, 0.03, 0.01, and an ex- 
ample cluster with N — 10 4 particles and A = 0.1 is 
shown in figure ^ At regular intervals, the growth of 
the clusters was suspended and 10 5 probe particles were 
"fired at" the cluster: these probe particles were allowed 
to diffuse freely, one at a time, until they contacted the 
cluster, at which point their location was recorded and 
the particle was deleted. In this way the growing prop- 
erties of the clusters (such as the penetration depth and 
multipole moments) were estimated. The code used is a 
direct descendent of that of Meakin [l4|, which in turn 
builds on the computational tricks of Ball & Brady j2^| 
to speed up computation. Whilst the code is truly lattice 
free, the smallest step size that a particle was allowed to 
take was set equal to one particle radius. Comparisons in 
2 dimensions between this and an algorithm which uses a 
much smaller minimum step size have shown that any ef- 
fect that this has on the results is the same order as the 
noise in the measurements attributable to inter-cluster 
variability |24| . 



I. FRACTAL DIMENSION AND CORRECTION 
TO SCALING EXPONENTS 

We calculate the effective value of fractal dimension 
from the local slope of the average radius of deposition 
vs number of particles N according to 



D = 



ln(iY 2 )-ln(iVi) 



ln(R dep (N 2 )) -la(Rdep(Ni)) 



(2) 



where the properties are measured at two different cluster 
sizes Ni and N 2 - For the clusters grown, the value of the 
fractal dimension is shown in figure [21 The dimensions 
estimated for each value of A appear to be converging to 
a common value of D = 2.50 ± 0.01 which is consistent 
with previous computational estimates |13( . The data for 
A = 0.03 and A — 0.01 are less well converged, and the 
results could be made more accurate with data for larger 
N. 

We now consider the correction to scaling exponents 
for DLA in three dimensions. In two dimensions it has 
been shown that no property of the cluster has a slower 
correction to scaling than the relative penetration depth 
and suggested that all properties should show influence of 
this slowest correction Q • To measure the slowest correc- 
tion to scaling exponent of a quantity, we used differential 
plots which proved effective for DLA in two dimensions 
12]. For some converging quantity Qn, which displays a 
single correction to scaling (eq. , then 



dQ 



N 



d\n(N) 



-v(Q N - Q c 



(3) 



so a plot of d fi^y) against Q should exhibit a straight line 
with slope —v, intercepting the x-axis at the asymptotic 



value Qco- The differential is approximated by 
d\n{N) ' 

and the error in the differential by 



ln (f) 



(4) 



dQ \ ^vHQn^ + ctHQn,) 



dhi(N) J 



In 



(ft) 



(5) 



where u{Qn) is the standard error in Qn- 

To characterise the DLA clusters we measured the 
multipole powers, since the corresponding multipole mo- 
ments provide an orthogonal set which may be used to 
totally describe the growing properties of the clusters. In 
three dimensions the multipole moments are estimated 
by (see |2^|, chapter 4) 



— 1,71 



l,m \ 



(6) 



using n probe particles which contact the cluster at 
(ri,9i,(fii) for i = 1 to n. We normalised the multipole 
power as 



Pi 



Em=-1 \ll,m\ 



(2l + l)Rf s 

where the effective radius is in turn given by 
1 



R. 



eff 



(7) 



(8) 



The definition of R e Q is such that it gives the radius of 
spherical target of equivalent absorption strength to the 
cluster. Note that for each measurement we used the 
centre of charge as origin, meaning zero dipole moments 
and hence P\ — 0; otherwise there is confusion between 
cluster shape and drift of the cluster centre (albeit the 
latter is rather negligible). 

The first and important feature of our differential 
plots is the indication of limiting asymptotic values 
Qoo (corresponding to d -^ N ^ — 0) which are universal, 
independent of the level of noise reduction. This is shown 
for the relative penetration depth in fig. El multipole pow- 
ers, Pi - P5 being shown in fig. El and the relative variabil- 
ity of extremal cluster radius in fig. [S] The universality 
of the asymptotic values of all these plots is strong in- 
dication that the limiting distribution of cluster shape is 
universal. 

The slopes of these same plots indicate the correction 
to scaling exponents, which also appear to be universal 
with respect to the noise reduction. Figure shows the 
measured values of the correction to scaling exponents for 
each of the quantities plotted (and all multipole moments 
measured). The exponent for the quadrupole power P 2 
is significantly different from the exponents for — and 
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P% . There is no quantity which shows a slower correction 
to scaling than -g*— , which suggests that the result found 

by Somfai el al. pfl in some sense also applies to clusters 
grown in three dimensions. 

The values of the correction to scaling exponents are 
least precise for the highest multipole moments, as these 
are most sensitive to the fine structure of the cluster. In- 
triguingly, the asymptotic value of the relative penetra- 
tion depth for DLA clusters in three dimensions is equal 
within measurement error to that of the two dimensional 
case: see figure[3]here and figure 3 in [lq. The correction 
to scaling exponent for is around | the value of the 
exponent for clusters grown in two dimensions, indicating 
that DLA clusters in three dimensions are considerably 
slower to mature. 



cluster with the same radius as a non-noise reduced clus- 
ter of N particles. Our data below show that this is a 
systematic underestimate, so that each value of A corre- 
sponds to a more severe noise reduction than expected. 
Figure shows the two point correlation function for 
DLA clusters grown at different noise reductions. The 
graphs have been shifted vertically so that all the curves 
collapse to a single line. From the shift factors used, we 
estimate the effective noise reductions to be 



.4 

A = 
A = 



= 1 
0.3 
0.1 



eff 



= 1 
0.19 
0.05 
e eff = 0.012 
A = 0.01 e eff = 0.0034. 



e 
£ eff 

.eff 



(13) 



A = 0.03 



II. FIXED POINT NOISE REDUCTION 

It is clear from differential plots such as figure |3] that 
the noise reduction "controls" the slowest correction to 
scaling. For low values of A this correction to scaling is 
strongly reduced, and we may write the behaviour of this 
correction as follows 



Qn = Qoo + C(A)N- 



(9) 



Other corrections to scaling need not depend on A, but it 
is quite evident that the amplitude of the slowest correc- 
tion to scaling is strongly affected by it. If there is some 
value of A for which C(A) is zero, then this value of the 
noise reduction would corres pon d to the fixed point of a 
renormalisation scheme (see |l0|). The plots for Pi — P± 
suggest that the fixed point noise reduction is A? < 0.01 
and plots of -J*— , -P5 are unclear as to the value of the 
fixed point noise reduction. Hence, one estimates that 



A f < 0.01. 



(10) 



The noise reduction of the fixed point can also be es- 
timated from the asymptotic properties of DLA clusters 
using the Barker & Ball ^3 formula 



D 



2 I oR ex t 



\ Re 



(11) 



where R ex t is the extremal cluster radius (the radius of 
the furthest cluster particle from the seed particle) and 
SR ex t is the cluster to cluster variability of R ex t- From 
figure [S] the asymptotic value of the relative variability 



of extremal cluster radius is 
This leads to 



8 Re 



0.0064 ±0.0016, 



0.032 ± 0.004. 



(12) 



which is close to the estimated value of A' . 

For a noise reduction of A, one would naively assume 
that it would require of order N / A particles to grow a 



Hence one concludes that the fixed point noise reduction 
in equation ITU1 should be adjusted to 



e f < 0.0035. 



(14) 



The values for ef and e* differ by a factor of 2, demon- 
strating that the identification process is open to some 
errors. If, as indicated by the results in figure [5] a sin- 
gle slowest correction to scaling exponent does not con- 
trol the scaling of all parameters, then a renormalisation 
scheme which is based on a single parameter (the noise 
reduction) may be inaccurate and agreement between e* 
and e* is not expected to be perfect. 



III. CONCLUSION 

The growth off-lattice of noise-reduced diffusion- 
limited aggregates in three dimensions has been consid- 
ered and shown to exhibit universality with respect to 
noise reduction. The fractal dimension is found to be 
D = 2.50 ± 0.01 which agrees with previous computa- 
tional and mean field j26J estimates. The penetration 
depth scales with the radius, and the asymptotic value 
of the relative penetration depth is -p^— = 0.122 ± 0.002 
which overlaps the value found for clusters grown in two 
dimensions |l2j. The convergence of the multipole pow- 
ers provides a very strong indication that DLA cluster 
growth, in three dimensions and off-lattice, converges to 
a universal distribution of cluster shapes. 

The relative penetration depth exhibited the slowest 
correction to scaling, jV°- 22±0 03 , and some multipole 
powers and also the relative fluctuations in extremal 
radius exhibited correction to scaling exponents which 
could be consistent with the same. However not all quan- 
tities exhibit the influence of the slowest correction to 
scaling and in particular the convergence to scaling of 
the dipole power, P2, is significantly faster than that of 
either or P 3 . 

Reducing the input noise by factors up to 100, by 
growing clusters in shallow bumps, clearly reduces the 
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amplitude of the leading correction to scaling. This sup- 
ports in three dimensions the idea of Barker & Ball |10| 
that the intrinsic fluctuation level is the physical origin 
of that slowest correction to scaling. We estimated the 
fixed point noise reduction to be e* ~ 0.0035 and this 
is close to the value estimated using the Barker & Ball 
|10| formula in terms of relative fluctuation in extremal 
radius. 

Taken together our results support the hypothesis that 
isotropic DLA in three dimensions approaches a sim- 
ple ensemble of statistically self-similar clusters, with a 
rather slow approach to scaling which is associated with 



the level of local geometric fluctuation. From this point 
of view, a quantitative model of the convergence of that 
fluctuation level appears to be the outstanding challenge 
in understanding isotropic DLA (in any dimension). An- 
other important challenge for three dimensions, for which 
work is in progress, is the role which material anisotropy 
can play |27|. 
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FIG. 1: A DLA cluster grown in three dimensions with N = 
10 4 particles and noise reduction factor A = 0.1, where the 
different shading indicates a different time of deposition on 
the cluster. 
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FIG. 2: The measured fractal dimension of DLA clusters, 
estimated by taking the local slope of Rd ep - The dimension 
appears to be converging to a dimension of D = 2.50 ± 0.01, 
universal with respect to value of noise reduction A. 
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FIG. 3: Differential plot of Relative Penetration Depth — 
against its own derivative with respect to \nN. The intercept 
with zero derivative indicates the asymptotic value of 

for infinite N is given by I = 0.122 ± 0.002 and the 

common limiting slope of the plots indicates a correction to 
scaling exponent of v = 0.22 ±0.03. Growth at different levels 
of noise reduction A > 0.03 is consistent with universal values 
of asymptote and exponent, whilst A — 0.01 appears to start 
and remain close to the 'fixed point' value of R ^ . 
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FIG. 4: Differential plots for the first four (non-trivial) multipole moments Pi — P$. All of the plots exhibit universal asymptotic 
values, corresponding to the extrapolation to zero derivative, and universal limiting slope corresponding to their correction to 
scaling exponent. The correction to scaling exponents were estimated by eye as: v{P?) = 0.32 ± 0.02, v(Pa) = 0.24 ± 0.03, 
v{Pa) = 0.26 ± 0.06, and f(Pb) = 0.29 ± 0.05, where the errors represent the maximum believable error. 




FIG. 5: Differential plot of the relative fluctuation in ex- 
tremal radius, Ms^. The asymptotic value is S ^"^ | = 

0.032 ± 0.004 which leads to an estimate of the fixed point 
noise reduction of e* = 0.0064 ± 0.0016. 




FIG. 6: The correction to scaling exponents obtained from 
differential plots of different multipole powers, and the two 
other quantities marked. The exponent for the dipole power 
P2 is significantly different from the exponents measured for 
~rf^ and P 3 . 



9 




L 



FIG. 7: The two point correlation function for clusters grown 
in three dimensions, scaled by effective noise reduction fac- 
tors, e ul . This noise reduction is chosen so that a data col- 
lapse is seen for small L. 



